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Abstract: The class of Chaplygin gas models regarded as a candidate of dark energy can be 
realized by a scalar field, which could drive the variation of the fine structure constant a during 
the cosmic time. This phenomenon has been observed for almost ten years ago from the quasar 
absorption spectra and attracted many attentions. In this paper, we reconstruct the class of Chap- 
lygin gas models to a kind of scalar fields and confront the resulting Aa/a with the observational 
constraints. We found that if the present observational value of the equation of state of the dark 
energy was not exactly equal to —1, various parameters of the class of Chaplygin gas models are 
allowed to satisfy the observational constraints, as well as the equivalence principle is also respected. 
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1. Introduction 



The possibility of varying the fundamental constants over cosmological time-scale has been studied 
for many years |jj . Among these fundamental constants, the time- variation of the fine-structure 
constant a is deserved to study both from the experimental and theoretical point of view. The 
observational evidence of time- varying a is firstly from the quasar absorption spectra and reported 
in 2001 by Webb et al., and now there are some other independently results can be used to constrain 
the variation of a, see [[TJ], Q and ||. The variation of a can be due to many possible reasons, one 
of them is that there is a scalar field coupled to the gauge field. In 1982, Bekenstein first introduced 
the exponential form for the coupling, which in practice can be taken in the linear form coupling 
between a scalar field and the electromagnetic field to explain the variation. 

As we mentioned above, there are several observational constraints on the variation of a as the 
function of redshift z, namely 

Aa{z) _ a(z) - a ^ 
a oiq 



where ao = a(0) denotes its present measured value. After Webb et al.'s report, the Oklo natural 
fission reactor found the variation of a with the level |Aa/a| < 10 -7 at the redshift z = 0.14 H], 
1 4]. The computation of 187 Re half-life in meteorites gives |Aa/a| < 10 -7 at the redshift z = 0.45 
|2|. The absorption line spectra of distance quasars suggests Aa/a — (—0.543 ± 0.116) x 10~ 5 
between z — 0.2 and z — 3.7 [||. And the recent detailed analysis of high quality quasar spectra 
gives the lower variation Aa/a — (—0.06 ± 0.06) x 10~ 5 over the redshift z = 0.4 — 2.3 and 
Aa/a = (0.15 ±0.43) x 10" 5 over the redshift z = 1.59-2.92 ||. In the following, we only consider 
the upper limit for the variation of a, so we have chosen the conservative constraint | Aa/a\ < 10~ 6 
over the redshift z = 0.4 — 3. The limit from the power spectrum of anisotropies in the Cosmic 
Microwave Background (CMB) is |Aa/a| < 10~ 2 at z = 10 3 Q. At last, the most ancient data 
from Big Bang Nucleosynthesis is |Aa/a| < 10 -2 over the redshift z = 10 8 - 10 10 §, §). 

Quite independently, observations like Type la supernovae, CMB and SDSS et al. have strongly 
confirmed that our universe is accelerated expanding recently caused by an unknown energy com- 
ponent called dark energy. Experiments have indicated there are mainly about 73% dark energy 
and 27% matter components in the recent universe, but so far people still do not understand what 
is dark energy from fundamental theory. The best candidate seems the cosmological constant in- 
cluding the vacuum energy, but it suffers the fine-tuning and coincidence problems. In order to 
alleviate these problems, a kind of scalar models called quintessence is needed to explain the origin 
of the dark energy. Thus it is natural to consider that quintessence or other type of scalar field 
models could be responsible for the time variation of a and it is a possible way to to distinguish 
dynamical dark energy models from the a cosmological constant. For the recent progress on the 
varying alpha, see || and references therein. 

Another candidate class of dark energy model is called the Chaplygin gas model Jl(J , which has 
been developed for many years, see (ll[] , [[[2) and Jl3| . Since Chaplygin gas can be always described 
by a scalar field with a effective potential, which is a process called reconstructing a scalar field from 
Chaplygin gas. So, it is natural to consider the time variation of a driven by Chaplygin gas which 
realized by a scalar field. In this paper, we first briefly review how a scalar field (quintessence) to 
drive a time-varying a in the next section. In Section 3. various Chaplygin gas models are reviewed 
and we also reconstruct them to scalar fields. In Section 4. we show the variation of a driven by 
Chaplygin gas and the conclusion is given in the last section. 
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2. Varying alpha from quintessence 



Let us consider the following action 

S = \j ^x^f~g TZ + J d A x^} L^{4>) + J d 4 x^9 C F ((f>, F^u) + S m , (2.1) 

where is the Lagrangian density for the quintessence field (fr minimally coupled to gravity as the 
following 

C+^l&W^-Vtt), (2.2) 

and S m is the action for ordinary pressureless matter. Here, C F is the Lagrangian density for an 
electromagnetic field F^ u coupled to quintessence field as 

C F {<t>) = -\B F {<j>)F^F»" , (2.3) 

where B F {(fr) allows for the evolution in <fr and Bp {(fro) = 1, where the subscript represents the 
present value of the quantity. The effective fine structure constant depends on the value of (fr as 

-W (2 ' 4) 



and thus we have 

Aa _ a — an 1 — Bp{(fr) 
a a B F {(fr) 



(2.5) 



We will consider a spatially flat Friedmann- Robertson- Walker (FRW) universe with the metric ds 2 = 
dt 2 — a 2 {t)dx 2 and assume the scalar field is homogeneous during the evolution of universe, then 
the energy density and pressure of the scalar field are 

p* = Y + v ^> p* = y- v W> (2 - 6) 

where a dot denotes the derivative with respect to the cosmic time t. And the equations of motion 
are 

H = -\(p m + ^j, (2.7) 
p m = -3H Pm , (2.8) 
^.SH^ffl, (2.9) 



subject to the Friedmann constraint 



H 2 = l( Pm + P* ) , (2.10) 



and the solution to eq.( p.8[ ) is simply p m = p m oa~ 3 . In fact, the equation of motion (2.S) for <fr 
should be added a term proportional to F^F^ and the derivative of Bp. However, such a term can 
be safely neglected due to the following reasons. First, the derivative of Bp actually corresponds to 
the time derivative of a, which is very small when we consider the equivalence principle constraints. 
Second, the statistical average of the term F^F^ over a current state of the universe is zero. 
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3. C haply gin gas models 



3.1 Chaplygin gas 

There exist an interesting class of dark energy models involving a fluid known as a Chaplygin gas 
@, which can explain the acceleration of the universe at later times and its equation of state is 



A 

P = , 

P 



(3.1) 



which can be obtained from the Nambu-Goto action for d-branes moving in a (d + 2)-dimcnsional 
spacetime in the light-cone parametrization. With the equation of state (3.1) the energy conserva- 
tion law d(pa 3 ) — —pd(a 3 ) can be integrated to give 



A + 4" , 

a 6 



(3.2) 



where B is an integration constant. By choosing a positive value for B, we can find that p ~ yB/a 3 
when a is small (a <C {B/A) l / & ) and p ~ —p ~ V~A when a is large (a 3> (B/A) 1 / 6 ). Thus, at earlier 
times when a is small, the gas behaves like a dust (pressureless) and it behaves as a cosmological 
constant at late times, thus leading to an accelerated expansion. In a generic situation, there is an 
intermediate phase, in which it looks like a mixture of a cosmological constant with a "stiff" matter 
= P)- 

One can obtain a homogeneous scalar field 4>(t) with its potential V(tf>) and Lagrangian density 
(2.2) to describe the Chaplygin cosmology by setting the energy density and pressure of the field 
(2.6) equal to that of the Chaplygin gas and we find 



B 



aS^A + B/a 6 



V{4>) = - y/A + B/ a e 



A 



y/A + B/aP 



(3.3) 
(3.4) 



By using the Friedmann equation ( 2.10| ), we get the variation of tf> in terms of the integration of a 
from eq.( 3.15| ): 

W ' l<! (3-5) 



oo = 



x a(Aa 6 + B) 1/4 (3H 2 ^mo + (Aa e + B)^/ 2 ) 1 / 2 ' 



where tt m o = /9 m o/(3i?Q) is the energy density parameter of matter and we have set ao — 1. The 
equation of state of Chaplygin gas is w = p/ p = —A/(A + B/a 6 ) > — 1, then we get 



B = — 



A(l + w ) 



Wq 



(3.6) 



and from the Friedmann equation (2.10), we obtain 



3H 2 (i-n m0 ) = VaTb. 



(3.7) 



From eq.(3.6) and (3.7) we get 



.4 



-w 



3^0 (1 — ^mo) 



B = (l + w Q ) 



3^(1 - 



nd) ) 



(3.8) 
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and then eq.(3.5) becomes 



0-0o = V3(l + w ) [° (1 + w - woa 6 )- 1 ^ (-^- + (1+ w -woa 6 ) 1 / A ^ - 



(70-(70-l)a 6 ) 1/4 



To 



1/4 



tan 



\ 



1/4 / 2 

To V r +V 



tanh 



1/4 / ; j 

To v r +v 



V 



1/2 

r - To 



r + 7 ' 



/ 



(70-(70-l)a 6 ) 1/4 



y=i 



(3.9) 



where 70 = 1 + w Q , r — f2 m o/(l - ^mo) < 1 and y = (70 - (70 - l)a 6 ) 1/4 . Therefore, if w = -1, 
is a constant during the evolution of the universe as the cosmological constant. When 

1/6 



a > 



1 + w — r A 
w 



we can neglect r in eq.(3.5) and get 



— = tanh 



1/2 
vTo , 



(70-(70-l)a 6 ) 1/4 

y=i 



(3.10) 



(3-11) 



3.2 Generalized Chaplygin gas model 

Although Chaplygin gas provides an interesting possibility for the unification of dark matter and 
dark energy. However, it have to face some problems to explain some current observations such as it 
leads to the loss of power in CMB anisotropics. This problem could be alleviated in the generalized 
Chaplygin gas model proposed in ref. Jll| (also see |Q) with equation of state 



A 



(3.12) 



where < a < 1 and when a = 1 it reduce to the pure Chaplygin gas ( |3.lD . Together with energy 
conservation law d(pa 3 ) — —pd(a 3 ), it gives 



A 



D 



3(1+<t) 



(3.13) 



where B is an integration constant. Hence, we can see that, for small and large a it behaves like a 
dust and a cosmological constant respectively, but in the intermediate phase, it looks like a mixture 
of a cosmological constant with a "soft" matter whose equation of state is p — op which can be 
obtained by expanding the pressure (3.12) and energy density (3.13) in subleading order: 



p ~ -A 1 *" 



oBA i+" 



,-3(X+ff) 



p ~ A^+" + 



BA~ — 



,-3(i+<0 



(3.14) 



1 + cr ' 1 + cr 

One can also reconstruct a minimally coupled scalar field to mimic the behavior of generalized 
Chaplygin gas by identify its energy density ( 3.15 ) and pressure ( |3.12j ) to that of the scalar field 

L> — , (3.15) 



,3(1+*)/ ^ + 



V(<f>) = 



2^ a 3(l+<r) + B 



2a 3 (!+ CT ) A 



(3.16) 
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By using the Friedmann equation ( 2.10|) , we get the variation of (j> in terms of the integration of i 
from eq.(3.15|): 



'SB da 



(3.17) 



The equation of state of generalized Chaplygin gas is w = p/p — —A/(A + B /a 3 ( 1+a )) > —1, then 
we get 

A(l+w ) 



B = -- 



w 



which is the same as eq.(3.£) and from the Friedmann equation (2.10), we obtain 



3iT 2 (l - O m0 ) = {A + B)-> 



From eq.(3.18) and (3.19) we get 



A = ~w 



3-^0 (1 — ^mo) 



1 + (T 



B = (1 + w ) 



1 + <T 



(3.18) 



(3.19) 



(3.20) 



and then eq.(3.17) becomes 



0o = ^3(1 + wo) J ( 1 + w - w a 3(1+ff) 



'2(1 + <t) 



r + ( 1 + w - w a 3(1+<T) 



-1/2 



da 
a 



'To 
3 Ji 



To 



_( 70 _l )Q 3(i+-)l2TT+^T 



2y 1+a dy 



( y 2 + r )l/2 ( y 2(l+a) _ 7() J ' 



(3.21) 



where 70 = 1 + W , r = n m0 /(l - Q m0 ) and y = [ 7o - (70 - l)a 3 < 1+<T )] 2<1+ "» . Therefore, eq.(|33l|) 
can be analytically calculated when r could be neglected, namely: 



a > 



7o - r 
7o-l 



3(l+<r) 



and then we get 



(3.22) 



2V ~ 3 .tanh-U 



3(1 + 0-) 



1/2 
7o 



3.3 Modified generalized Chaplygin gas model 



(3.23) 



Another candidate for the generalization of the Chaplygin gas called modified generalized Chaplygin 
gas or modified Chaplygin gas model p3| is characterized by the following equation of state 



p = Ap 



B 



(3.24) 



where A, B and a are constants and < cr < 1. Thus, it looks like a mixture of two kinds of fluids, 
one with equation of state p = Ap and the other one being the generalized Chaplygin gas. From 
eq.( |3.24| ), one can see that it reduces to generalized Chaplygin gas when A = and to the perfect 
fluid if B = 0. Again, together with energy conservation law, it gives 



where C is an integration constant. Then, for small scale factor a, it behaves like a dust (if A = 0) 
or radiation (if A = 1/3) with equation of sate p = Ap and energy density p = C T +° r a _3 ' 1+ " 4 ', wile 
for large a, it behaves like a cosmological constant. In the intermediate phase, it corresponds to the 
mixture of a cosmological constant and a kind of fluid with equation of state p — (a + A + Aa)p, 
which can be obtained by expanding the pressure (3.24) and energy density ( 3.25| ) in subleading 
order: 



B \~ 



1+A 
B \~ 



1 + A 



a + A + Aa 
B 



B 



-TTS Ca ~3(l+A)(l+a) 



1 + A 



1 + A / 

TTo 7 C , a -3(1+A)(l+ CT ) 



1 + CT 



1 + CT 



(3.26) 
(3.27) 



We also reconstruct a minimally coupled scalar field to mimic the behavior of generalized Chaplygin 
gas by identify its energy density ( 3.25 ) and pressure ( |3.24 ) to that of the scalar field 



(1 + A)C 



v(4>) 



n 3(l+A)(l+a) / B , C 

25(1 + 4)-i a 3(i+A)(l+a) + (1 - A)C 



2a 3 ( 1 + A )( 1 + CT ) 



c 



1+A ~ a 3 ( 1 + A )( 1 +") 



(3.28) 



(3.29) 



By using the Friedmann equation ( 2.10| ) , we get the variation of in terms of the integration of a 
from eq.(HJf): 



£>o = 



y/3(l + A)C da 



The equation of state of generalized Chaplygin gas is 

p „ 5(1 + A) 



"'"p ^ B + (1 + A)Ca- 3 ( l +W+°) 



> -1 



where we have used eq.q3.24j) and ( |3.25j ). From eq.(|3.31|), we get 

B(l+w ) 



and 



, dw 



°~ {l + A)(A-w ) 
3BC(l + A) 3 (l + cr) 



5=0 [s+(i+A)cr 



3(l + (T)(l + too)(A-tco) 



1/2 ' 

(3.30) 
(3.31) 

(3.32) 
(3.33) 



where we have used eq.( 3.32 ) and z = a 1 — 1 is the redshift. From the Friedmann equation ( 2.1C ), 
we obtain 



3^0 (1 — ^mo) — 



B 



1 + A 



C 



Using eq.([3~3|), |3~33|) and ( jO^ ), we get 



B = (A - 7b) 



3B (1 — O m0 ) 



1 + w 

A 



1 1 + fT 



3B~o(l — n m o) 



(3.34) 



(3.35) 
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where we have defined 70 = 1 + wq and 



A= l + A 



70 • 



3(1 + (7)70 

Noticed that the case A = corresponds to w' = 370(1 — 70) (1 + c)- Then eq.( 3.30D becomes 



(3.36) 



0o = V^To 



.4 



7o + (A- 7o )a 3A ( 1 +-) 



To 
3 h 



^ + (A-^)a^(l + „)] 2IT+57 



2y 1+ °dy 



da 
a 



(y 2 + ra 3 ( A - 1 A V2 (Ay 2 ^) - 7o ) 



(3.37) 



where r = Q m o/(l - O to0 ) and y 
cally calculated when r could be neglected, namely 



70 + (^-70> 3A(1+g) 



Therefore, eq.(3.37) can be analyti 



r < y 2 a 3(ji 1} , or a > 



7o 



a 3(i-l)(l+ff) _ ^ r i+o 



70 -A 



(3.38) 



and then we get 



- ^tanh-if^ M 

3(l + cr)Vi \ V 




3A(1 + <t) ] 2(l + cr) 



(3.39) 



4. Varying alpha from Chaplygin gas models 

Although the form of the coupling between a scalar field and the electromagnetic field can be very 
complicated |l4| ] , in general, the observational results have indicated that the variation of a is small 
and Bp can be approximated as a linear form in practice. Therefore, in this paper, we will take 
such a approximation 

S F (0)=l-C(0-0o), (4.1) 

which corresponds to the choice of e = r = and q = 1 for the parameters in ref. |Q. From the 
tests of the equivalence principle, the coupling is constrainted to be |£| < 10~ 3 . Then, the variation 
of a is given by 



Aa 



|C(0-0o)|. 



(4.2) 



b f {4>) 

Here, we tried some different values of £ to make all the constraints that mentioned in the introduc- 
tion section be satisfied. The variation of a is presented in Fig.|l|, Fig.|| and Fig.|| for the Chaplygin 
gas model, the generalized Chaplygin gas model and the modified generalized Chaplygin gas model 
respectively. 

From Fig.|l|, one can see that all the constraints are respected for ( < 1.1 x 10~ 6 in the case 
of the Chaplygin gas model. For the case of generalized Chaplygin gas and a given value of £, 
the variation of a is getting smaller and smaller when the parameter a becomes small in the same 
situation. In other words, the smaller a is, the larger upper bound of £ is, see Fig]|. Finally, in 
the case of modified generalized Chaplygin gas, the variation of a becomes large when the present 
running of the equation of state w' Q is large. Thus, the upper bound of £ should be smaller to satisfy 
constraints, see Fig.||. 
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Figure 1: Chaplygin gas: The variation of a during the evolution of the universe is plotted, namely, 
log |Aa/a| vs. -log(l + z). Here we have used wo — 0.99, Q m o = 0.27 and the solid, dashed and dot-dashed 
cures correspond to £ = 1.1 x 10~ 6 , £ = 5.0 x 10~ 6 and £ = 0.2 x 10~ 6 respectively. Only the curves that 
not overlaps the gray areas are phenomenologically viable. 




-Log(l+z) -Log(l+z) 



Figure 2: Generalized Chaplygin gas: The variation of ot during the evolution of the universe is plotted, 
namely, log |Aa/a| vs. -log(l + z). Here we have used wo = 0.99, fi m o = 0.27, £ = 1.1 X 10~ 6 and the solid, 
dashed and dot-dashed cures correspond to a = 1.0, a — 0.5 and a = 0.1 respectively. Only the curves that 
not overlaps the gray areas are phenomenologically viable. 



5. Conclusions 



In this paper, we have reconstructed the class of Chaplygin gas models to a kind of scalar field and 
study the variation of the fine structure constant a driven by it. This phenomenon was found since 
ten years ago and attracted many attentions. The resulting |Aa/a| as a function of the redshif z 
is presented in Fig.|l], Fig.|] and Fig.||. We only consider the case of linear coupling between the 
scalar field and the electromagnetic field, because the variation of a is much small. The results 
indicate that if the present observational value of the equation of state of the dark energy was not 
exactly equal to —1, various parameters of the class of Chaplygin gas models are allowed to satisfy 
the observational constraints, as well as the equivalence principle is also respected since it requires 
the constant £ is much smaller than 10 -3 in all the case. 

For the generalized Chaplygin gas, there is a parameter < a < 1 in eq.(3.12). We find that 
when the smaller a is, the larger upper bound of ( is. In the case of modified generalized Chaplygin 
gas, the upper bound of £ becomes more restricted when the running of equation of state w' is 
large. It is worth further studying, since the variation of fundamental constants during the cosmic 
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-10 -8 -6 -4 -2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

-Log(l+z) -Log(l+z) 

Figure 3: Modified generalized Chaplygin gas: The variation of a during the evolution of the universe is 
plotted, namely, log|Aa/a| vs. -log(l + z). Here we have used wo = 0.99, Q m o = 0.27, £ = 0.2 x 10 -6 , 
a — 1.0 and the solid, dashed and dot-dashed cures correspond to A = (or w' = 0.06), w' — 0.12 and 
Wq ~ 0.24 respectively. Only the curves that not overlaps the gray areas are phenomenologically viable. 



time is a very interesting area. 
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